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We evaluate the electron-electron interaction corrections to the electronic thermal conductivity 
in a disordered conductor in the diffusive regime. We use a diagrammatic many-body method 
analogous to that of Altshuler and Aronov for the electrical conductivity. We derive results in one, 
two and three dimensions for both the singlet and triplet channels, and in all cases find that the 
Wiedemann- Franz law is violated. 
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I. INTRODUCTION 

The effect of electron-electron interaction on the elec- 
trical conductivity of disordered systems has been exten- 
sively investigated both theoretically and experimentally 
over the past two decadeslii. There are two main types of 
correction to the Drude electrical conductivity which can 
have similar magnitudes and temperature dependences: 
weak localization and interaction effects. Weak localiza- 
tion is due to interference between pairs of time-reversed 
scattering trajectories of electrons from impurities; this 
effect can occur in a non-interacting system. Interaction 
effects are due to the increased effective electron-electron 
interaction strength due to incomplete screening by dif- 
fusively moving electrons. Experimentally the two effects 
can be distinguished by applying a magnetic field; weak 
localization is suppressed whilst interaction effects are 
not. 

In comparison there has been relatively little work 
done on the thermal conductivity, either theoretically 
or experimentally. This is largely because thermal con- 
ductivity is hard to measure in low-dimcnsional sys- 
tems at low temperature, and it is difficult to separate 
the electronic and lattice contributions. It is therefore 
doubly hard to observe the disorder-driven corrections 
to the electronic thermal conductivity. Moreover there 
have been theoretical predictions'^''' that the Wiedemann- 
Franz law holds, which allows one to deduce the thermal 
conductivity, k, directly from the electrical conductivity, 

where ks is Boltzmann's constant, e is the electronic 
charge, T is the temperature, and Lq is known as the 
Lorenz number. Clearly it is important to establish 
whether the Wiedemann-Franz law is valid in interacting 
disordered systems; if it is, there is nothing to be gained 
from measuring the thermal conductivity in addition to 
the electrical conductivity. 

In this paper we evaluate the interaction corrections 
to the thermal conductivity in one, two and three dimen- 
sions in both the singlet and triplet channels. We find 



that the Wiedemann-Franz law is violated in all cases. 
Our results are presented in detail in Table 1. We have 
not included a result for the singlet channel in one di- 
mension since this is just one third of the correspond- 
ing triplet channel result with F replaced by an effec- 
tive singlet channel interaction Fs = KgO^ In (D/ig/T), 
as explained in the figure caption. Each term consists 
of two pieces: the first piece leads to violation of the 
Wiedemann-Franz law; the second piece does not. The 
exception is the triplet channel in two dimensions which 
does not possess a violating piece. Note that in each case 
both terms have the same parametric dependences; it is 
only the constant prefactors which are different. 

Our calculation is the exact analogue of the original 
Altshuler- Aronov calculation for electrical conductivity^. 
There have been several previous attempts at this 
calculation"*'^'^; however these are in disagreement with 
each other, and we believe all of them to be incorrect. 
The fact that it has taken so long (over twenty years) 
to get the correct result for the thermal conductivity is 
due to three inherent difficulties in the problem. First, 
the heat-current operator is not uniquely defined, and 
some choices of definition will be renormalized by the 
electron-electron interaction^. Second, the heat-current 
operator has opposite electron-hole parity to the charge- 
current operator, which leads to the presence of an extra 
"heat-drag" diagram for the thermal conductivity which 
is vanishingly small for the electrical conductivity. (We 
use the term "heat-drag" by analogy to the Coulomb 
drag effect between two layers of electrons; this effect is 
described by the same diagrarr^). Third, spurious ultra- 
violet divergences occur in the diagrammatic approach to 
thermal conductivity; these can be understood and eval- 
uated within the framework of divergent series theorjiiS 
(their origin is due to illegal series expansions in the 
derivation of the perturbation theory) . 

Before we proceed to the details of our calculation, we 
present a short history of the field. We will first con- 
sider effects which do not require interaction, such as 
weak localization and the Anderson transition. Chester 
and Thellung'^ used an exact eigenstates approach to 
show that the Wiedemann-Franz law should hold in a 
non-interacting disordered system, independent of the 
strength of impurity scattering. Strinati and Castellaniii 
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TABLE I: Interaction corrections to the electronic thermal conductivity in one, two and three dimensions in both the singlet 
and triplet channels. In the above equations T is temperature, D is the diffusion constant, F is the effective interaction in the 
triplet channel, and K2 is the inverse screening length in two dimensions. Note that h — 1 and fcs = f in these calculations. 
We do not include a result for the singlet channel in one dimension as it is just one third of that for the triplet channel with F 
replaced by the effective interaction parameter for the singlet channel, Fs ~ t^^a? In (Dk^/T), where K3 is the inverse screening 
length in three dimensions, and a is the wire width. The inverse screening lengths in two and three dimensions are given by 
k| — 8'KN{0)e'^ and K2 = iTTN{0)e^ , where A'^(O) is the single-spin electronic density of states in the appropriate dimension. 



used a Ward identity construction to argue that the 
Wiedemann-Franz law holds all the way to the Ander- 
son transition. Kearney and Butcher^^ used the exact 
eigenstates approach of Chester and Thellung to deduce 
the weak localization correction to thermal conductivity. 
This effect was later seen in the experiment of Bayot et 
aU^, which measured the electrical and thermal magneto- 
conductance, which were found to obey the Wiedemann- 
Franz law. Enderby and Barnesi^ later pointed out that 
the Wiedemann-Franz law is not obeyed at the Anderson 
transition - the previous treatments all use the Sommer- 
feld expansion, which is not valid close to the Anderson 
transition. However the only effect is to decrease the 
constant in the Lorenz number from 7r^/3 « 3.29 to 2.17 
for a conductivity exponent 1/ = 1. The overall con- 
clusion for the non-interacting effects is that there is no 
extra information in the thermal conductivity that is not 
present in the electrical conductivity. In the weak local- 
ization regime the Wiedemann-Franz law holds; close to 
the Anderson transition it is only modified by a simple 
numerical factor, so all critical properties are correctly 
predicted by it. 

It would therefore seem that the only possibility of in- 
teresting behavior lies with the interaction corrections 
to thermal conductivity. Castellani et al"* predicted that 
the Wiedemann-Franz law would hold even for the in- 
teracting disordered system. These authors evaluated 
the dynamic energy-energy correlation function in the 
interacting system using a skeleton graph analysis, and 
derived the thermal conductivity from this. Since their 
method of calculation is somewhat different from ours, 
it is not obvious to us why they obtain a different re- 
sult. Livanov et alS later directly calculated the inter- 
action corrections using a quantum kinetic equation ap- 



proach, and predicted that the Wiedemann-Franz law is 
violated in all dimensions. They appear to have evalu- 
ated some interaction contributions correctly, but to have 
missed other contributions of the same order of magni- 
tude; for example, in two dimensions, they predict a log- 
arithmic increase in k/T at low temperature rather than 
the correct logarithmic decrease. We believe that the 
results obtained initially from their diagrammatic for- 
malism are correct, and that they have then erroneously 
thrown away some terms believing them to be smaller in 
magnitude. Arfi'' performed the equivalent calculation 
in the Matsubara formalism and again found violation 
of the Wiedemann-Franz law. However this calculation 
has errors relating to all three of the difficulties referred 
to previously: (a) the Matsubara form of heat current is 
used, but the heat current is then erroneously renormal- 
ized by interaction; (b) the "heat drag" or Aslamazov- 
Larkin-like diagram is omitted; (c) the final results are 
parametrically larger than those for electrical conductiv- 
ity by factors of l/(rr)^, where r is the elastic scatter- 
ing time, which appears to be due to incorrect treatment 
of spuriously divergent terms, as will be described later. 
The aim of the present work is therefore to resolve the dis- 
crepancies between previous calculations and to present 
the first consistent calculation of the interaction correc- 
tion to thermal conductivity. 



The remainder of the paper is organised as follows: 
in section 2 we derive the Drude contribution to elec- 
tronic thermal conductivity; in section 3 we present an 
outline of the calculation of the interaction corrections to 
thermal conductivity; finally in section 4 we analyse our 
results and draw conclusions. 
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where N{0) is the single-spin electronic density of states 
at the Fermi surface. We see that the response function 
is isotropic, so we drop the spatial indices. We can then 
perform the e; sum to obtain 



k,e 



FIG. 1: Feynman diagram for the Drude contribution to the 
heat-current response function. The momentum and Mat- 
subara frequency dependence of the two electron lines are 
labelled. 



Qhh = -16^3^(0)i^r3^(/ + l/2-p/2)2 
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= -—N{0)DT^p^-p), 



(7) 



and we finally extract k using Eq. (|2Jl to get the Drude 
result 



Ko = —N{0)DT = — - — 



(8) 



II. DRUDE THERMAL CONDUCTIVITY AND 
THE WIEDEMANN-FRANZ LAW 

Before we proceed to the calculation of the interaction 
corrections, we demonstrate the simplicity and elegance 
of the Matsubara approach to evaluating thermal conduc- 
tivity by deriving the Drude result. The thermal conduc- 
tivity is obtained from the imaginary time heat response 
kernel, Qhh{i^p), by analytic continuation from positive 
Bose Matsubara frequencies, fip — 2nTp, 
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The Feynman diagram for the Drude thermal conductiv- 
ity is shown in Fig. 1. The solid lines are disordered 
electron Green functions 



G{k,iei) 



1 



iei - 6 + 27Sgn(e;) 



(3) 



where e/ — 2ttT{1 + 1/2) is a Fermi Matsubara frequency, 
= /c^/2m — /i is the electronic excitation spectrum, and 
r is the elastic scattering time. The black dots represent 
heat-current vertices, which are given by 



(4) 



in the Matsubara frequency representation that we use. 
The heat current kernel is then given by 

<(i^^p)-2r5:5:^^[i(e. + f7,/2) 

ei k 

G{k,iei)G{k,iei+i%) (5) 

Performing the k-integral, we only obtain a non-zero 
result if the frequencies e/ -I- Op and e; have opposite 
sign, which means that the e; is restricted to the range 
— Jlp < £; < 0. Upon changing the sign of ti we obtain 
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(ei-np/2)\ (6) 



where n is the electron number density and we have 
used the Einstein relation 2N(0)D = nr/m. From the 
corresponding Drude formula for electrical conductivity, 
(To = ne'^T/m, we see that the Wiedemann-Franz law is 
obeyed (note that h = 1 and fc^ = 1 in our calculations). 

This diagrammatic technique offers a simple proof that 
the Wiedemann-Franz law is obeyed for a non-interacting 
disordered system with arbitrary disorder strength. For 
any diagrammatic contribution to electrical conductiv- 
ity, fj, there is a corresponding contribution to thermal 
conductivity, k. Moreover the corresponding expressions 
differ only in the form of the current vertices, yielding a 
factor kakfse^/m? for tr, and kakf3[i{ei + Q.p/2)Y /m? for 
K. The only essential difference between k and a then 
lies in the frequency sums; the ratio of these two sums 
is independent of disorder and leads to the Wiedemann- 
Franz ratio. In particular, the weak localization correc- 
tion to thermal conductivity is obtained directly from the 
Wiedemann-Franz law. Note that this proof relies on the 
Sommerfeld expansion since we are linearising our energy 
integrals about the Fermi surface. Interaction effects can 
violate the Wiedemann-Franz law since the presence of 
an interaction line can alter the Matsubara frequencies 
at the two current vertices. 



III. THE INTERACTION CORRECTIONS 

In this section we calculate the interaction corrections 
arising from the singlet and triplet interaction channels— 
- we do not evaluate Cooperonic contributions since these 
are expected to be small for a system with repulsive inter- 
actions. In the singlet channel the dominant contribution 
arises from small energy and momentum transfers be- 
tween electrons. This is dominated by the bare Coulomb 
interaction, which takes the form 
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takes the form 



2NiO)Vtiq,iUn) = F 



Dq^ 



{F + l)Dq^ 



(13) 



where F = 2N{0)V^ and is the bare interaction in the 
triplet channeli^. Note that F includes Fermi liquid cor- 
rections in the absence of disorder (it is only unscreened 
with respect to the disorder). F may be determined from 
the measured paramagnetic spin susceptibility in the ex- 
perimental system of interest 





(E) 



(F) 



FIG. 2: Feynman diagrams for the interaction correction to 
the heat-current response function. Diagrams A-E are struc- 
turally identical to the diagrams which contribute to elec- 
trical conductivity. Diagram F is the additional "heat-drag" 
diagram which must be included to obtain a consistent result. 



where d is the dimensionality, and a is a measure of the 
transverse width in a quasi-one dimensional wire. The 
disorder-screened singlet interaction then takes the form 
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Voiq)~^+niq,iuJn)' 
where the polarization operator n(g,ia;„) is given by 

Dq^ 
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Iliq,iLU„) = 2N{0) 



Dq^ 



(11) 



If any integrals we obtain involving Vs{q, ii-Un) are conver- 
gent at small momentum, q, we can ignore the Vo{q)^^ 
term in Eq. (|10l) in comparison to Il{q,iu!n)', Vs{q,iWn) 
then takes the universal form. 
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In the triplet channel the dominant contribution arises 
from momentum transfers of the order of 2kp, and the 
unsceened triplet interaction can be treated as a con- 
stant. The disorder screened triplet interaction then 
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1 + F 
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If we first calculate the contribution from the singlet 
channel, we can then obtain the corresponding result for 
the triplet channel by replacing V, by Vt and multiplying 
by 3 - the extra factor of 3 arises because there are 3 
times as many modes in the spin-one channel than in the 
spin-zero channel. 

Applying standard perturbation theory we find that 
the leading order interaction corrections are given by the 
Feynman diagrams shown in Fig. 2. These diagrams are 
identical to those considered by Altshuler and Aronov 
for the electrical conductivitji^, with the exception of di- 
agram (F) which we call the "heat-drag" term. The cor- 
responding diagram for electrical conductivity is a factor 
{T/EfY smaller than the other diagrams due to cancella- 
tion of electron and hole charge currents; the heat current 
has opposite electron-hole parity and thus electron and 
hole heat currents reinforce each other. Consequently di- 
agram (F) is of the same size as the other diagrams in the 
case of thermal conductivity. The use of the Matsubara 
heat current vertex defined in Eq. ensures that no 
interaction renormalization of vertices is needed, greatly 
simplifying the calculation with respect to other choices 
of heat current operator (obviously the final results must 
be independent of this choice). 

Details of the calculation of the contributions from 
diagrams A-F to the heat-current response function, 
Qhh{i^), are presented in the appendix (note that in 
the following we suppress the subscripts on Matsubara 
frequencies Jlp — 2'KTp and LOn — 2'KTn for convenience) . 
It is shown there that Qhh{0) vanishes, as it must for 
internal consistency, and that kT may be written as the 
limit of the expression 
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as n tends to zero. As explained in the appendix, this 
statement is to be interpeted in the sense of first contin- 
uing iD,p to real frequencies, and then taking the limit 
of real frequency going to zero. The second line of Eq. 
(|15|l is just the interaction contribution to the electrical 
conductivity multiplied by the factor 7r^T^/3e^. In other 
words, this is the contribution to kT predicted by the 
Wiedemann-Franz law. The other two terms thus lead 
to violation of the Wiedemann-Franz law if they yield a 
non-zero result (as we indeed find that they do). The ex- 
pressions in Eq. H15|l may now be evaluated in one, two 
and three dimensions for both the singlet and triplet in- 
teractions to yield the results listed in Table 1. Note that 
all three terms in Eq. H15|l give results having the same 
parametric form but with different constant prefactors; 
we must therefore evaluate all of them to get a correct 
final result. This statement is not quite true in two di- 
mensions because of the presence of logarithmic terms, 
and we should keep the most singular terms in this case. 

IV. CONCLUSIONS 

We have calculated the interaction corrections to ther- 
mal conductivity in the diffusive regime of a disordered 
conductor. Our main result is that the Wiedemann-Franz 
law is violated in all dimensions; the predicted interac- 
tion results have the expected parametric dependences, 
but different numerical coefhcients. For example, in the 
singlet channel in two dimensions, the actual logarithmic 
correction is half that predicted by the Wiedemann-Franz 
law. Unfortunately the experimental work of Bayot et 
aU4 is the only work we know of on thermal conduc- 
tivity in disordered conductors, and this only isolates a 
weak localization contribution. We hope that further ex- 
periments will be performed in this area, and that the 
interaction correction be observed as well as the weak lo- 



calization correction. In particular the two-dimensional 
system would seem to be a promising one to investigate. 
Bayot et al's^'* disordered graphitic system showed weak 
localization effects whose magnitude was roughly 20% of 
the Drude term in the electronic thermal conductivity at 
a temperature of 2.9K. If it were possible to cleanly ex- 
tract the phonon term, one could look at the electronic 
term as a function of both temperature and magnetic 
field, and distinguish weak localization and interaction 
effects. This would then experimentally settle the ques- 
tion of whether the Wiedemann-Franz law is violated in 
an interacting disordered system. Another sensible quan- 
tity to investigate experimentally would be the thermal 
Hall conductivity, which arises solely from electrons (al- 
though there can be phonon drag effects). In future work 
we intend to derive the interaction corrections to this 
quantity. 
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APPENDIX: DIAGRAMMATIC 
CONTRIBUTIONS 

In this appendix we list the complete set of diagram- 
matic contributions for archival purposes. This should 
allow interested researchers to reproduce our results in 
detail. The total contribution to the heat-current re- 
sponse function, Qhh{i^), from all the diagrams in Fig. 
2 is given by 
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In the above expression we have hsted the contribu- 
tions according to which diagram they are obtained from; 
the additional numerical label following the letter refer 
to different sign configurations of Matsubara frequencies 
which are possible within the same diagram. For exam- 
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pie, in diagram (A), there are three possible sign configu- 
rations whose contributions we denote by (^41), (^2) and 
{A3) . To proceed further we next perform the sums over 
the Fermi Matsubara frequency, e. After some simplifi- 
cation we then obtain the following result for Qhh{i^)- 



V{q, 



-4N{0)Dd^pTYJ2' 

u)>Q q 

+16N{0)DT Y Y^'^('^ - ^) 

4 



m{0)Ddc,pT Y Y^^^ 



Dqaqf3V{q,uj) 



0<u}<Q q 



{Dq^ + ujf 



ui>Q. q 



{Dq^ + Oj)[{Dq^ + uj)^ - 



-NiO)DT yy\n' + Uco'n - An^T'n] , ^^"^f^^f 
3 UOq^ +uj)^(Dq^ +uj + n) 



+ -N{0)DT y y Uw' + div'n + 3u,n' - Att'T'u] , , ^^"^f 



u:>0 q 



{Dq^+ujf{Dq^+tj + nf 
I 



(ABC) 

(DE) 

(DE) 

(DE) 

(F) 



(A.2) 



Note that the contributions from diagrams A, B and C 
do not cancel each other as they do in the corresponding 
electrical conductivity calculation. At this point we can 
check that we have chosen a consistent set of diagrams 
to evaluate by setting the external frequency, f2 = 0, 
and checking that Qhh{^) = 0. This is a very powerful 
test that should not be omitted - it is very dangerous to 



merely calculate [Qhh{i^) — Qhh{^)]/^, especially when a 
large number of diagrams are involved. One always runs 
the risk of missing important physical processes; in fact 
we only became aware of the presence of the heat-drag 
term of diagram (F) when this check failed in its absence. 
Setting = here yields 
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where in the last step we have used the divergence theo- its bounding surface at infinity. Note that this check 
rem to convert the g-integral into a surface integral with 
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not only gives us confidence that all relevant diagrams 
have been included, but also that each diagram has been 
given the correct combinatorial factors. In the above 
derivation we have treated VQ{q)~^ as though it had no 
q-dependence. This is justified as the terms ignored are 
smaller in magnitude, and would be cancelled by higher 
order diagrams. In addition to the above proof that 
Qhh{0) = 0, we have also directly evaluated Qhh{0) from 
Eq. (jA.3|l in one, two, and three dimensions, for both 
singlet and triplet potentials, and found it to equal to 
zero. 

We can then find kT as the limit of Qhh{^p)/^p as 
rip tends to zero. Now of course this statement does not 
make sense since Hp is discrete - we should first analyti- 



cally continue to real frequencies and then take the limit. 
However it turns out to be possible to avoid performing 
this analytic continuation, and instead to continue ma- 
nipulating discrete sums. We can then extract the term 
that is proportional to Qp and discard terms proportional 
to higher powers of Qp. This approach is legitimate pro- 
vided that any operation we perform on f2p would carry 
over unchanged to the same operation on the correspond- 
ing real frequency after analytic continuation. The ad- 
vantage of this method is simply convenience in calcula- 
tion - no illegal operations occur, as we have checked by 
performing the analytic continuation first and then tak- 
ing the zero-frequency limit. This procedure leads to the 
expression for kT as the zero-frequency limit of 
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We can then perform the g-integral using the standard 
identification 
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and the formulae for the V{q, iuj) given in Eq. H1(J|I and 
Eq. (|13|l . After performing this g-integral we end up with 
terms which are infinite sums over powers of w. Provided 
that the sums are ultraviolet divergent, we can identify 
them as zeta functions via 

T ^ w*^' = T{2'r:Tf ^ = r(27rT)'=C(-fc)- (A.6) 



Any sum that is infrared divergent must be cut off cor- 
rectly since such a divergence is physical - such a situa- 
tion occurs in two-dimensions where we obtain a logarith- 
mic sum which is cut off at w ~ T at low frequency and 
~ 1/r at high frequency. We have used this particular 
divergent series trick as it allows a very direct evaluation 
of results. As a check of its legality we have recalculated 
the various terms using standard analytic continuation 
methods and obtained the same results - albeit after a 
lot more algebraic manipulation. 
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